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IRAS PSCZ VS. IRAS 1.2-JY MODEL VELOCITY FIELDS: A SPHERICAL HARMONICS COMPARISON. 
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ABSTRACT 

We have used the two IRAS redshift surveys, 1.2-Jy (Fisher et al. 1995) and PSCz (Saunders et al. 2000), to 
model the linear velocity fields within a redshift of 8 000 km s" 1 and have compared them in redshift space. The 
two velocity fields only differ significantly in their monopole components. The monopole discrepancy cannot be 
solely ascribed to shot-noise errors and incomplete sky coverage. The mismatch seems to arise from incomplete- 
ness of the PSCz catalog at fluxes < 1.2 Jy. The 1.2-Jy and PSCz higher order velocity multipoles, particularly 
the dipole and quadrupole components, appear to be consistent, suggesting that the dipole residuals found by 
Davis, Nusser and Willick (1996) when comparing 1.2-Jy and Marklll velocity fields probably originates from 
the Mark III velocities calibration procedure rather than from uncertainties in the model velocity field. Our results 
illustrate the efficiency of the spherical harmonics decomposition techniques in detecting possible differences 
between real and model velocity fields. Similar analyses shall prove to be very useful in the future to test the 
reliability of next generation model velocity fields derived from new redshift catalogs like 2dFGRS (Colless et 
al. 2001), SDSS (York et al. 2000), 6dF and 2MRS. 

Subject headings: cosmology — theory - galaxies large-scale structure of the Universe — large-scale dynamics. 



1. INTRODUCTION 

In the framework of linear gravitational instability theory and 
linear biasing a comparison between the predicted and observed 
peculiar velocities allows us to measure the so called (3 param- 
eter 

0=-^-. (1) 
b 

i.e. a combination of the mass density parameter, Q m , and the 
linear bias parameter, b, which relates the mass overdensity 
field S to the density fluctuations in galaxy counts 6 g , through 
S g = bS. Several different techniques have been used to compare 
measured peculiar velocities (taken from all available catalogs: 
Marklll Willick 1997a, SFI Haynes 1999, ENEAR da Costa et 
al. 2000, SEcat Zaroubi 2002) to predict velocities that, in most 
cases, have been obtained from two redshift surveys of IRAS 
galaxies: the 1.2-Jy and the PSCz catalogs. 

Most analyses have investigated the consistency between 
model and observed velocities, while only a few of them have 
been devoted in comparing different datasets or in studying the 
consistency of model velocity fields. All recent works show 
a general consistency between models and data (Willick et al. 
1997, hereafter WSDK; Willick and Strauss 1998; Sigad et 
al. 1998; Dekel et al. 1999; Branchini et al. 1999, hereafter 
B99; Nusser et al. 2000; Zaroubi et al. 2002) with one no- 
ticeable exception represented by the work of Davis, Nusser 
& Willick (1996, hereafter DNW). Their analysis showed that 
the velocity residuals between the Marklll peculiar veloci- 
ties and the 1 .2-Jy predictions display a significant spatial cor- 
relation. Subsequent analyses by WSDK and Willick and 



Strauss (1998) were performed using the same dataset, a sim- 
ilar model velocity field but a different comparison technique, 
called VELMOD, which allows for independent calibrations of 
Marklll velocities in each of the Marklll sub-catalogs. As 
a result the observed and model velocity fields were brought 
into better agreement, at least within a redshift distance s = 
7 500 km s" 1 and returned a calibration inconsistent with the 
original one performed by Willick et al. (1997a). 

While these results suggest that the DNW mismatch orig- 
inates from calibration problems in the Marklll catalog, the 
possibility that systematic errors in the 1.2-Jy model velocities 
also contribute to the discrepancy has been completely over- 
looked. In particular, we ignored whether the velocity fields 
predicted from the deeper IRAS PSCz catalog provides a better 
match to the Mark III velocities than the shallower 1 .2-Jy cata- 
log, considered by DNW. 

The main goal of this work is to check whether the discrep- 
ancies between the 1.2-Jy and Mark III velocities can be as- 
cribed, at least to some extent, to model uncertainties. For 
this purpose we compare both the overdensity and the radial 
velocity fields modeled from the redshift space distribution of 
1.2-Jy and PSCz galaxies. This work is meant to stress the 
usefulness of spherical harmonics decomposition techniques in 
assessing the adequacy of a model velocity (or density) field 
through velocity-velocity (or density-density) comparisons. 

In § 2 we describe how to compare velocity and overden- 
sity fields inferred from different redshift surveys. The galaxy 
redshift catalogs used in our analyses are described in § 3. In 
sections 4 and 5 we analyse the two model overdensity fields 
and radial velocity fields. Our main results are discussed in § 6 
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and our main conclusions are presented in § 7. 

2. SPHERICAL HARMONIC COEFFICIENTS DECOMPOSITION 

The relations between an arbitrary scalar field ^>(s) and the 
real-valued spherical harmonic coefficients ipi m (s) are 



oo / 



(2) 



(3) 



1=0 m=-l 

where ^; m (s){Z = 0, .., oo; \m\ < 1} denote the well known real- 
valued spherical harmonics, defined as in Baker et al. (1999, 
hereafter BDSLS) and Bunn (1995) (see also Jackson 1999). 
For a given / these functions are normalized to the value of 34o 
at the North Galactic Pole. 

The amplitude of the Z-multipole ipi(s) is defined by the sum 
in quadrature over the (21+ l)-ip[ m 's, where {m = -l,...,l}: 
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m=l 
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(4) 



2.1 . Peculiar velocity and overdensity fields from the 
distribution of galaxies in redshift space 

Nusser & Davis (1994) show that in the linear regime the 
peculiar velocity field is irrotational not only in real space but 
also in redshift space and thus can be expressed as the gradient 
of a scalar velocity potential: v(s) = -V$(s). At a given redshift 
s, <3?(s) and 5(s) can be expanded in spherical harmonics and 
related to each other through a modified Poisson equation: 

1 ( 2 v 1 l(l+l)$ lm /3 A ldlntf. 

(5) 

where <j)(s) and <5/ m (s) denote the sample selection function 
and the spherical components of the overdensity field, respec- 
tively. Here prime expresses d/ds. Solving this differential 
equation requires computing the redshift space density field 
from the observed galaxy distribution on a spherical grid us- 
ing Gaussian cells of approximately equal solid angle. These 
cells are equally distributed in longitude (64 bins) and are cen- 
tered at the Gaussian-Legendre 32 point quadrature formula in 
the range -1 < cos(fr) < 1, where b is the Galactic latitude. 
In this work we measure redshifts in the Local Group frame, 
and use 52 Gaussian radial redshift bins out to a distance of 
s = 1 8 000 km s" 1 . When comparing two or more galaxy cata- 
logs the size of the radial bin is set by requiring equal spatial 
resolution and the minimal shot noise constant throughout the 
volume. We do this setting the distance between the centers 
of the Gaussian cells at distance s equal to the average galaxy- 
galaxy separation, a(s) = [n</>(s)] -1 / 3 , in the sparser catalog. The 
mean number density is estimated as 
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1 



(6) 



where the sum is over the galaxies contained within the volume 
V, which, in this work, is a sphere of radius 8 000 km s" 1 . 
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FIG. 1. — The 1.2-Jy smoothing width as a function of redshift. The dots 
represent the value of a„ at centre of the 52 redshift bins. 



Hernquist and Katz (1989) pointed out that there are two pos- 
sible ways to estimate the smoothed density field: the "scatter" 
and "gather" approches. In the former approch each particle is 
distributed in space, and the density estimate at a given point is 
the superposition of individual smoothing spheres. In the lat- 
ter case one defines a smoothing lenght at a given point and 
weights the particle in its neighborhood by the resulting kernel. 
The Gaussian-smoothed galaxy overdensity field at a grid cell 
centred around s„ is given by: 



l+S(s n ) = 
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(7) 



where the sum is over the N galaxies in the catalog, which 
clearly is a "gather" approach. The smoothing width at a red- 
shift s, used in this paper is <r n = max[320, a(s)] km s" 1 . Since 
in this work we compare the PSCz catalog with the sparser 
1.2-Jy catalog, we use the <r„ obtained from the 1.2-Jy catalog 
shown in Fig. (1). This smoothing scheme is tailored to keep the 
shot-noise uncertainty roughly constant throughout the sample 
volume and is similar to the optimal Wiener filtering procedure 
used by Lahav et al. (1994) and Fisher et al. (1995b). More- 
over, the our kernel has a compact support, which simplifies 
the computational task. 

Eqn. (5) is valid in the linear regime, where a one-to-one 
mapping between distance and redshift is guaranteed. This as- 
sumption is not valid in high density regions, such as cluster of 
galaxies, where shell crossing may occur triple valued-regions 
may apear. To deal with the latter effect, we adopt the same 
procedure as Yahil et al. (1991) and collapse the fingers of God 
associated with the six richest clusters in the sample listed in 
Table 2 of that work. 

In this work we will consider galaxies within s max = 20000 
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kms" 1 with the spherical harmonics decomposition limited to 
'max =15. Outlying s max galaxies are assumed to be distributed 
uniformly according to the mean number density of the parent 
catalogs. 

3. DATASETS 

The IRAS PSCz catalog contains rj 15 500 IRAS PSC galax- 
ies with a 60 fim flux larger than 0.6 Jy. The average depth of 
this survey is ks 100 h~ l Mpc . We restrict our analysis to a sub- 
sample of 1 1 206 galaxies within 20 000 km s" 1 from the Local 
Group (LG, hereafter). The areas not covered by the survey 
amount to 16.0% of the sky and are preferentially located near 
the galactic plane, in the so-called Zone of Avoidance (Saun- 
ders et al. 2000). The empty areas have been filled-in with the 
cloning procedure described in Branchini et al. (1999) . 

The 1.2- Jy catalog (Fisher et al. 1995a) contains 5 331 IRAS 
galaxies with a 60 /im flux larger than 1.2-Jy within 20000 
kms" 1 . This catalog has a slightly larger sky coverage of (« 
87.6%) and a smaller median redshift (« 84 h~ l Mpc). The 
same B99 filling technique has been used to restore all sky cov- 
erage. 

Fisher et al. (1995b) have found that uncertainties in filling - 
in the empty areas of the 1 .2-Jy survey do not appreciably affect 
the spherical harmonics coefficients with I < 15. Similar con- 
clusions are also valid for the IRAS PSCz survey (Teodoro et 
al. 1999). 

The selection functions of the two surveys, that quantify the 
probability of a galaxy at distance s, have been determined as 
maximum-likelihood fits to the smooth function 

0(s)=j (tO (iS*) ,s>H r s , (g) 
[ 1 , otherwise, 

where r = s/Hq 1 is expressed in /r'Mpc (Yahil et al. 1991). 
The best-fit parameters, determined using the galaxies within 
8 000 kms -1 and assuming a flat universe with Vt„, = 1, are 
shown in Tab. (1). The quantity r s in Eqn. (8) accounts for 
the incompleteness of faint galaxies in the innermost regions of 
the catalogs (Rowan-Robinson et al. 1990, Yahil et al. 1991). 

4. THE IRAS PSCz AND 1.2- Jy FIELDS OVERDENSITY 

The first step of our analysis consists of applying Eqn. (7) to 
compute the IRAS 1 .2-Jy and PSCz overdensities onto a spheri- 



cal grid of radius 18 000 km s 
ing scheme given by Eqn. (7). 



4.1. 



and l m ., x = 15 with the smooth - 



Overdensity Maps 

In Figs. (2—4) we show the Aitoff projections of the PSCz (top 
panels) and 1 .2-Jy (middle panels) galaxy overdensity fields in 
three spherical shells centred at s = 1 000, 3 000 and 5 000 km s" 1 , 
respectively. The thickness of each shell is defined by the Gaus- 
sian width at the distance of the shell. All maps are plotted in 
Galactic coordinates. The bottom panels show the map of the 

'Throughout this article we write Hubble's constant as Hq = lOOA 
kms -1 Mpc -1 




FIG. 2. — PSCz, (top), 1.2-Jy (middle) and residual (bottom) overdensity 
fields in the radial shell centred at s = 1 000 km s . The thick-line contour 
indicates 8 = or <5 re s = (bottom). The light- and dark-gray contours show 
over- and under-dense regions, respectively. In the top and middle panels the 
contours are equally-spaced by A<5 = 1.0. The bottom panel the density resid- 
uals are shown in steps of A<5 re s = 0.125. 



overdensity residuals S les (s) = Sq.(,(s)- 61.2(3), where the sub- 
scripts 0.6 and 12 refer to PSCz and 1.2-Jy surveys, respec- 
tively. In all plots the thick line indicates the zero overdensity 
(or residual) contour. The overdensity (residual) contour spac- 
ing AS (A(5, es ) depends on the shell's depth and is indicated in 
the Figure captions. Detailed cosmographic studies of the 
structures in the nearby universe have already been performed 
using both the distribution of 1.2-Jy (Strauss and Willick 1995) 
and PSCz galaxies 2 and will not be repeated here. We only 
stress that the main cosmic structures like the Virgo Cluster (/ « 
240°, Z?«75°, s= 1000 kms" 1 ), the local void spanning the re- 
gion 330° </ < 120°, -90° <b< 30° at s= 1000 kms" 1 , the 
Hydra-Centaurus complex (/« 270°, /?« 20°, s = 3000 kms" 1 ) 
and (Z « 300°, b « 20°, s= 3000 kms" 1 ) and the Perseus- 
Pisces supercluster (I « 140°, b « -25°, s = 3000 kms" 1 ) are 
detected in both surveys. 

The residuals, calculated at the same gridpoints as the den- 

2 See http://www-astro.physics.ox.ac.uk/~wjs/pscz.html 
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Table 1 

Selection function parameters for IRAS catalogs. 



Catalog /i™„, Ngal' d a (3 1 r s b r* n 

[Jy] [A _1 Mpc] [A _1 Mpc] [/z 3 Mpc" 3 ] 



PSCz 


0.6 


13 364 


0.54 


1.80 ... 


6.0 


87.00 


6.12 


•10- 


-2 


PSCZ 0.7 


0.7 


11 170 


0.57 


1.80 ... 


6.0 


76.35 


6.38 


•10" 


-2 


PSCZ 1.2 


1.2 


5519 


0.43 


1.86 ... 


5.0 


50.40 


4.62 


•10" 


-2 


1.2-Jy 


1.2 


6010 


0.43 


1.86 ... 


5.0 


50.60 


4.67 


•10" 


-2 


-PiS Cz evo. 


0.6 


13 364 


0.99 


3.45 1.93 


6.0 


76.02 


5.76 


•10" 


-2 



a Includes synthetic objects. 

b Parameter kept constant in the likelihood maximization. 

c The selection function parameters have been taken from Canavezes et al. (1999) 
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Table 2 

The parameters of the 81.2-80.6 and K1.2-H0.6 linear regressions. 



Regression 


Ngrid AW 


A a 


B 




^1.2-^0.6 


112 552 156.36 


-0.011 ±0.014 


1.025 ±0.036 


0.78 


"1.2 -«0. 6 


112 552 < 156.36 


60.0±5.9 


1.040 ±0.019 


0.88 



Notes: 

° For the u — u comparison A is expressed in km s . 

Column 2: iVgrid, the number of gridpoints used for the regression; 

column 3: A'dof- the number of independent volumes; 

column 4: A, the zero-point offset in the linear regression, and its 1 - 

column 5: B, the slope of the linear regression, and its 1 — a error; 

column 6: Xeff/^dof from the linear regression. 



sity fields, seem to indicate that, as expected, the largest residu- 
als are located near the galactic plane where uncertainties in the 
filling-in procedure dominate the errors. In the inner shell, over- 
densities in the 1.2-Jy galaxy distribution appear to be larger 
than the PSCz ones while in the most external shell the situa- 
tion is reversed. 

4.2. Density Residuals 

A way of quantifying possible systematic differences be- 
tween the two fields is that of performing a point-to-point com- 
parison. In absence of systematic errors we expect the two 
IRAS overdensity fields to be linearly related through 

6ia=B6oj6+A (9) 

with A = and B = 1, 5q& and S1.2 are estimated at Cartesian 
gridpoint positions. A zero-point offset A 7^0 indicates discrep- 
ancies in the mean densities of the two catalogs while a slope 
B ^ 1 reveals the presence of systematic errors. If the values of 
60.6 and 61.2 are independent and their errors are Gaussian then 
we can compute the % 2 statistics 



v grid 



(S^-A-BSp.e) 2 



(10) 



'1.2 



where ero.6 an d C1.2 indicate the typical errors in 5 and N^d is 
the total number of points in the comparison. However, since 
the gridspacing (250 km s" 1 ) is smaller than the smoothing length 
(> 320 km s" 1 ), not all values of <$ .6 (or $1.2) are independent. 
Indeed, the effective number of independent points used in the 
X 2 statistics, A/ eff is given by 



^grid ^grid 

N£ = A^gnd" 2 £ £ expC-r*^.), 
<=i 7=1 



(11) 



where ry is the separation between gridpoints i and j and R s j 
is the smoothing radius of the Gaussian kernel at the gridpoint 



i. This expression is a generalization of Eqn. (12) in Dekel et 
al. (1993) to the case of a grid with variable mesh-size since 
R s i depends on the gridpoint position. The new statistics xls = 
(A/grid /A/eff) x 2 is approximately distributed as a x 2 with A/ dof = 
A4ff- 1 degrees of freedom and can be applied to infer the values 
of A, B and their errors. 

Fig. (5) shows the S1.2 vs. Sq.6 scatter diagram obtained us- 
ing all 112552 gridpoints with s < 8000 kms" 1 and \b\ > 8°. 
Within these limits the estimated overdensities are little affected 
by uncertainties in the filling-in procedure, and we expect that 



I to 1 > e° 

s < 8 000 kms"' 
A = -0.011 ± 0.014 (l-o): 
B = 1,035 1 0.036 (1-cr) 
a, ■■= Q. 16 




- 1 



FIG. 5. — 1.2-Jy vs. PSCz overdensity. Only 1883 randomly selected grid- 
points with \b\ > 8° and within s = 8 000 km s~' are shown in the scatter plot. 
The continuous line represents the linear best fit characterized by the zeropoint 
A and slope B indicated in the panel. 
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shot-noise dominates the total error budget. To estimate the 
shot-noise errors affecting both PSCz and 1.2-Jy overdensity 
fields we have generated 100 bootstrap realizations of the ob- 
served PSCz and 1.2-Jy surveys. These are obtained by replac- 
ing each galaxy (including the 'synthetic' galaxies used to fill- 
in the unobserved regions) with a number of objects drawn from 
a Poisson deviate with mean unity. Shot noise errors in 5 at the 
generic gridpoint were set equal to the dispersion over the 100 
realizations. The solid line in the plot is the best linear fit ob- 
tained from the minimization of x%f- As shown in Tab. (2), the 
slope of the line, B = 1.025 ±0.036, is consistent with unity, 
indicating no systematic mismatch. On the other hand the zero- 
point value A = -0.01 1 ± 0.014 shows no systematic difference 
in the mean density of the two IRAS catalogs. 

The scatter around the best-fitting line, as = 0. 16 is similar to 
average shot-noise error in the 1 .2-Jy density field (cr^ 2 = 0.15). 
The parameter S = x eff /A^ f = 0.78 indicates that bootstrap er- 
rors slightly overestimate the true uncertainties. 

4.3. Radial Density Contrast 

An alternative way to characterize possible discrepancies be- 
tween the two IRAS density fields is that of computing the dif- 
ferential radial density contrast in shells with identical thick- 
ness As for both galaxy distributions: 

p(s) = 1 X w(s, Si ,As) 
p Vnj^ cb(si) ' 1 

where the window associated to the spherical shell Ws p /p(s,Sj, 
As) is the Heaviside function H(s + As/2-Sj)H(sj + As/2-s), 
Sj is the redshift of the object i, V is the volume used to com- 
pute the mean number density n and the sum runs over the N p 
galaxies of the catalog. Fig. (6) shows the quantity p(s)/p for 
1.2-Jy (light, dotted histogram) and PSCz (thick, continuous 
histogram) galaxies within s = 15 000 kms -1 , computed in in 
spherical shells of thickness As = 200 km s" 1 . Although the 
general shapes of two radial density contrasts are similar (with 
a bump at s ~ 1 100 km s -1 due to Virgo and Fornax clusters, 
a dip at s ~ 3 500 km s" 1 induced by the Sculptor Void and a 
second prominent peak at s ~ 5000 kms" 1 , at the same red- 
shift as the Perseus-Pisces supercluster and the Great Attrac- 
tor), a discrepancy between the two curves exists in the range 
2 000 < s < 6 000 km s" 1 where the PSCz density is larger than 
the 1.2-Jy one. 

The two smoothed curves in Fig. (6) show the monopole 
terms of the two overdensity fields m= o as a function of 
radius, obtained from the spherical harmonics decomposition. 
Both curves interpolate the histograms rather well (hence show- 
ing the adequacy of our spherical harmonics analysis) and con- 
firm the existence of a PSCz vs. 1.2-Jy density mismatch be- 
tween the two catalogs. 

5. THE IRAS PSCz AND 1 .2-Jy PECULIAR VELOCITY FIELDS 

To perform a similar analysis on the 1 .2-Jy and PSCz model 
velocity fields we have integrated Eqn. (5) using (3 = 0.6, a value 



consistent with recent determinations based on several velocity- 
velocity comparisons (e.g. Nusser et al. 2000), a recent density- 
density analysis (Zaroubi et al. 2002) and the study of mean 
relative velocities of galaxy pairs (Juszkiewicz et al. 2000). 

5.1. Radial Peculiar Velocity Maps 

To obtain maps of the radial velocity field we have computed 
the gradient of the velocity potential projected along the line of 
sight: «(s) = s • V<J> at the same gridpoint positions used for the 
analysis of the density field. 

The top end mid panels of Figs. (7-9) are similar to those of 
figs. 2-A and show the radial velocity maps in the same shells. 
The bottom panels show the maps of the radial velocity residu- 
als M les (s) = UQ.e{s)- u\2(s). The thick line represents the zero 
radial velocity (residual) contour. The velocity (residual) con- 
tour spacing Au (A« res ) varies with the redshift and is indicated 
in the Figure captions. 

The main feature seen at all redshifts in both 1.2-Jy and 
PSCz radial velocity maps is the dipole pattern resulting from 
the LG reflex motion toward the CMB apex {I rj 276°, b « 30°, 
KogutefaZ.1993). 

Local deviations from this pattern arise from peculiar mo- 
tions, like the infall onto Perseus-Pisces superclusters at (s w 
3000 kms -1 , I « 120°, b w -40°). The largest velocity resid- 
uals are located near the galactic plane, but their contours are 
broader due to the intrinsic non-locality of the peculiar velocity 
field. The most striking feature is perhaps the fact that velocity 
residuals, that are positive across a large fraction of the sky in 
the first shell, becomes mostly negative in the remaining two 
shells, suggesting possible systematic errors in either PSCz or 
1.2-Jy velocity models. 




5000 10" 1.5x10" 

s [Kms- 1 ] 

FIG. 6. — Differential radial density contrast for PSCz (thick line histogram) 
and 1.2-Jy galaxies (thin line histogram). The monopole of the PSCz and 1.2- 
Jy density fields are shown with smooth solid and dashed lines, respectively. 
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FIG. 7. — PSCz (top), 1.2-Jy (middle) and residual (bottom) radial velocity 
fields in the spherical shell centred at s = 1 000 kms~' . The thick-line contour 
indicates Su = or A% e s = 0. The light- and dark-gray contours show regions 
characterized by outflows or infall motions. In the top and middle panels the 
contours are equally-spaced by Au = 200 km s _1 . The bottom panel shows 
the velocity residuals in steps of A<5 res = 25 km s _1 . 



5.2. Velocity Residuals 

To perform a more quantitative comparison between the two 
radial velocity fields we have repeated the analysis of § 4.2 and 
performed a linear regression between u j 2 and M0.6 by minimiz- 
ing 



A' 



grid / 



1.2" 



(«?.2 + B2 <0.6 



(13) 



where N gI a runs over the same points used for the density- 
density comparison. Linear velocities are obtained by integrat- 
ing the mass density fields over large scales causing them to 
be correlated on scales larger than the radius of the Gaussian 
kernel. This means that the number of independent gridpoint 
N e ff is smaller than the estimate given by Eqn. (11). Intrinsic 
large scale smoothing also affects the velocity error estimate, 
a u since uncertainties in modeling the galaxy distribution in the 
unobserved areas now add to shot noise errors. 




Same as in Fig. (7). The shell is now centred at s = 3 000 km s 1 and 
200 km s~ l and Ai( res = 50 km s _1 



FIG. 

the contour spacings are Au 



The shot noise errors have been computed using the same 
bootstrap resampling technique described in § 4.2. Uncertain- 
ties in the filling-in procedure have been evaluated from 20 in- 
dependent PSCz and 1 .2-Jy mock catalogs extracted from the 
iV-body simulations performed by Cole etal. (1997) simulating 
a volume of w 170/r'Mpc in an Einstein-de Sitter universe 
with a non-vanishing (Oa = 0.7) cosmological constant and 
density fluctuations normalized to the observed cluster abun- 
dance (Eke, Cole & Frenk 1996 ). Model velocity fields from 
each IRAS mock catalog were compared with velocities ob- 
tained from ideal, all-sky IRAS mock catalogs. The variance 
over 20 mocks at each gridpoint quantifies the uncertainties 
in the filling procedure. Total uncertainties were obtained by 
adding in quadrature the two errors. 

The results of the linear regression are summarized in Tab. (2) 
which shows the slope of the best fitting line B = 1 .040 ± 0.019 
and the value of the zero-point, A = 60.0 ±5.9 km s" 1 . While 
the slope is still consistent with unity, at the 2.1-cj confidence 
level, the zero point is significantly different from zero, hence 
corroborating the evidence for a mismatch in the average densi- 
ties of PSCz and 1.2-Jy catalogs. The dispersion around the fit 
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s = 5000 kms" 1 

FIG .9. — Same as in Fig. (7). The shell is now centred at s = 5 000 km s _1 and 
the contour spacings are An = 200 kms -1 and A% s = 50 kms -1 . 



is, a u = 88 kms" 1 is similar to the typical velocity error com- 
puted from the IRAS mocks and the bootstrap resampling cat- 
alogs (cr„ 06 = 71 kms" 1 and a Ul2 = 140kms _1 ). The parameter 
X eff /AW is less than unity which either indicates that errors 
are overestimated or that the value N e ff, taken from Eqn. (11), 
does significantly overestimate the true number of degrees of 
freedom. 

5.3. Velocity Multipoles 

The spherical harmonics decomposition allows us to bet- 
ter investigate and characterize possible discrepancies between 
1.2-Jy and PSCz velocity fields. Here we are only interested 
in the first three (I = 0, 1,2) components, corresponding to the 
monopole, dipole and quadrupole terms, which have already 
been investigated in previous analyses (e.g. BDSLS). It is worth 
stressing that the monopole and dipole components at a given 
redshift are only sensitive to the mass distribution within the 
corresponding distance and thus can be directly related to the 
structures within the volume considered, while the quadrupole 
components are also sensitive to the mass distribution beyond 
such a distance [see Eqn. (9) in Nusser & Davis 1994]. 
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FIG. 10. — 1.2-Jy vs. PSCz radial velocities computed at 1883 randomly 
selected gridpoint positions with \b\ > 8° and s < 8 000 km s _1 inferred radial 
velocity field (x-axis) versus that inferred from the 1.2-Jy survey (y-axis). The 
parameters of the best fitting line (A,B) the scatter a u are indicated. 



The monopole terms of PSCz (continuous line) and 1.2- 
Jy (dashed line) velocity fields are plotted in Fig. (1 1) as a func- 
tion of the redshift. The hatched and dashed regions represent 
the 1-(T uncertainty strip, computed from the mock catalogs and 
the bootstrap resampling procedure. 
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FIG . 11. — Monopole coefficients of 1 .2-Jy (dashed) and PSCz, (continuous) 
radial velocity fields. 1 — a errors are represented by the hatched (shaded) 
areas. 
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FIG. 12. — The velocity dipole coefficients, u\ m (s){m = —1,0, 1), inferred 
from PSCz (continuous) and 1.2-Jy (dashed) and their associated l-cr uncer- 
tainties. The three Galactic Cartesian components GX, GY, GZ are shown in 
the top-left, top-right and bottom-right panel, respectively. The total ampli- 
tude is shown in the bottom-left panel. The dotted line refers to the dipole of 
PSCz galaxies with flux /60^ m > 1-2 Jy- 



Both velocity monopoles display the same general features, 
as expected given the similarities between the density monopoles 
Sou(s) shown in Fig. (6). There is a significant radial in- 
flow around 5= 1 800 km s" 1 corresponding to the infall mo- 
tion toward the Virgo-Fornax complex, an outflow in the range 
1 800 < s < 4000 kms" 1 due to the presence of the local void 
and another infall beyond this radius, due to the Hercules, Hy- 
dra, Centaurus and Perseus-Pisces superclusters. Differences 
between the two velocity monopoles are detected with a signif- 
icance level larger than l-cr at very small radii (s < 1200km s" 1 ) 
and in the range 3500 km s" 1 < s < 7500 km s" 1 . 

Fig. (12) shows the total amplitude (lower left panel) and 
the three Cartesian components (remaining panels) of the two 
IRAS velocity dipoles. The total amplitude of the PSCz dipole 
(continuous line) is systematically smaller that the 1.2-Jy one 
(dashed line) but the discrepancy is well within the expected 
errors, in agreement with previous studies (e.g. Schmoldt et 
al. 1999). 

Finally, Fig. (13) shows the five components of the PSCz (con- 
tinuous line) and 1.2-Jy (dashed line) velocity quadrupoles, 
U2 m (s) (where m = -2 , . . . , 2) along with total amplitude (bottom- 
left panel). As for the dipole case, the two sets of quadrupole 
components agree to within l-cr apart from a small disagree- 
ment in the the m = +2 and m = -2 components beyond s = 7000 
kms" 1 . Indeed, beyond s w 6000km s" 1 the amplitude of the 
PSCz quadrupole decreases, while the 1.2-Jy one increases 
monotonically. 

It is worth noticing that our 1 .2-Jy velocity multipoles agree 




200 



200 



200 



7.5 2.5 

s/1000 [km s- 1 ] 

FIG. 13. — Quadrupole coefficients «2 m (s)(m = -2, ...,2) derived from 
PSCz (solid line) and 1.2-Jy (dashed) model velocity fields. Shaded regions in- 
dicate the I— a uncertainty regions. The dotted line represents the quadrupole 
components for PSCz galaxies with flux /go^,,, > 1.2 Jy. 



with those computed by BDSLS except that their multipole 
components exhibit sharper features at small radii {s < 500 
kms" 1 ). This mismatch derives from our use of a Gaussian 
smoothing kernel of radius 330 km s" 1 which is somewhat larger 
than that applied by BDSLS. 

6. DISCUSSION. 

Perhaps the most surprising result of our analysis is the dis- 
agreement between the PSCz and 1 .2-Jy velocity monopoles in 
the 1 800 < s < 4000 kms" 1 range. Before investigating the 
possible causes of this mismatch, it is worth noticing that none 
of our IRAS monopoles is consistent with an outward flow of 
magnitude w 400 kms" 1 at a distance of 7000 kms" 1 like the 
one detected by Zehavi et al. (1998) using a sample of 44 Type 
la supernovae. That result is also at variance with the velocity 
model obtained from the distribution of ORS galaxies obtained 
by BDSLS which means that, if confirmed, it would be difficult 
to reconcile with the Gravitational Instability scenario. 

A possible clue to understand the origin of the monopole 
mismatch is provided by the analysis of BDSLS that returns 
a velocity monopole somewhere between the PSCz and 1.2- 
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Jy ones in the range of interest. This seems to suggest that 
our monopole discrepancy originates from systematic errors in 
one of the IRAS catalogs, possibly some incompleteness of the 
largest catalog at faint fluxes. 

The results of Tadros et al. (1999) seem to confirm that 
PSCz catalog might be statistically incomplete at fluxes fwfjm < 
0.7 Jy. To test whether this is indeed the case we have extracted 
a brighter PSCz subsample (called PSCzo.i) by discarding all 
object with fu)^ m < 0.7 Jy and have repeated our spherical har- 
monics decomposition analysis. The results, shown in Fig. (14) 
indicate that the PSCzq.i velocity monopole (dotted line) basi- 
cally coincide with that of the complete sample. It is only when 
discarding galaxies with fluxes smaller than 1 .2 Jy that the ve- 
locity monopole (dashed line) agrees with the 1.2-Jy one. 

Another possibility of understanding the monopole mismatch 
is that of invoking a strong evolution of IRAS galaxies as a 
function of the redshift, which we have neglected in our calcu- 
lations so far. Roughly speaking evolution can be distinguish 
into two different types. Pure density evolution changes only 
the normalization of the luminosity function while its shape is 
invariant with redshift. Pure luminosity evolution on the other 
hand describes a possible change of the intrinsic luminosities of 
the sources but not their number density. Both forms of evolu- 
tion occur simultaneously in the Universe. For instance, galaxy 
merging reduces the number density of galaxies causing thus 
evolution. Besides, such events also cause different luminosi- 
ties of the final systems and hence leading to luminosity evo- 
lution. Other physical mechanism for galaxy evolution can be 
devised. Some evidences for a sizable evolution of IRAS galax- 
ies have been reported by a number of authors, although there 
seem to be much controversy about the amplitude of the effect 
(see Springel 1996 and references therein). 

Canavezes et al. (1998) expressed evolution in terms of a 
generalization of the selection function introduced by Yahil et 
al. (1991) given by 

^)e™i={ ^TijITF) ,s>H r s , (U) 

{ 1 , otherwise. 

Springel ( 1 996) showed that it indeed provides a rather accurate 
modeling of both IRAS data-sets. The estimated parameters us- 
ing a maximum-likelihood procedure are listed in Tab. (1) and 
are denoted by PSCz evo \- The effect of evolution are shown in 
Fig. (14) in which the velocity monopole term derived from 
the PSCz galaxy distribution in which evolutionary effects have 
been accounted for through the new selection function <fi(s) eYO i 
(dot-dashed line) is plotted against the unevolved model. The 
evolutionary effects are very minor, as was somewhat expected 
given the locality of our galaxy sample, and do not help in re- 
ducing the PSCz vs. 1.2-Jy velocity monopole mismatch. 

Overall, our analysis seems to indicate that the discrepancy 
between 1.2-Jy and PSCz velocity monopoles cannot be as- 
cribed to errors in the modeling of the selection function which, 
as we have seen, would not change the monopole terms appre- 
ciably. Catalog incompleteness at very low fluxes has to be 
ruled out too. We can only conclude that either the PSCz in- 
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FIG. 14. — Velocity monopoles for the 1.2-Jy (dark gray line) PSCz (light- 
gray), PSCzo.i (dot-dashed), PSCz 1.2 (long-dashed) PSCz evo j (dotted) sam- 
ples. 



completeness extend at objects with fluxes comparable to 1.2 
Jy or that the monopole mismatch does reflect a genuine differ- 
ence in the average density of faint vs. bright IRAS galaxies at 
s< 4000 kms"'. 

A second important result found in our analysis is that the 
1.2-Jy and PSCz dipoles agree within the l-cr errors. This 
means that the use of the PSCz rather than the 1.2-Jy catalog 
does not help in reducing the Mark HI-IRAS dipole residuals 
found by Davis, Nusser & Willick (1996). Indeed, both GX 
and GY components of the PSCz velocity dipole are much too 
large with respect to Mark III ones (see Fig 15 of Davis, Nusser 
& Willick 1996). Similar results were also obtained from the 
ORS vs. Mark III dipole velocity comparison (BDSLS). We 
conclude that the discrepancy between the Marklll and 1.2- 



Table 3 

Selection function parameters In the IRAS mock 

CATALOGS. 



Mock catalog 






r, 










[/r'Mpc] 


[/r'Mpc] 


PSCz 


0.53 


1.90 


10.90 


86.40 


1.2-Jy 


0.48 


1.79 


5.00 


50.40 
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Jy velocity fields cannot be alleviated by reducing the shot noise 
errors of the model velocity field. 

The agreement between the PSCz, 1.2-Jy and ORS dipoles 
is reassuring since it guarantees that uncertainties in model ve- 
locity fields arising from within the galaxy samples that derive 
from incompleteness in one (or more) catalogs, incorrect treat- 
ment of non-linear velocities or non-linear non-uniform biasing 
do not generate systematic errors apart from a the monopole 
mismatch in the range 1800kms -1 < s < 4000 kms" 1 . 

Unlike the case of monopole and dipole moments, possible 
differences between the 1.2-Jy and PSCz quadrupole compo- 
nents can also be attributed to incorrect determination and di- 
lute sampling of the density field beyond the sample. The fact 
that 1.2-Jy and PSCz velocity quadrupole agree within the ex- 
pected errors suggests that differences in the quadrupole com- 
ponents of the two fields can be understood in terms of shot 
noise errors and uncertainties in filling the empty regions. 

With this respect it is worth investigating whether the denser 
sampling of the PSCz catalog at large radii allows a better 
modeling of the large scale contribution to the model veloc- 
ity field. Indeed, when comparing the model 1.2-Jy velocity 
field to Mark III velocities in the framework of the VELMOD 
analysis, WSDK have found that velocity residuals exhibit a 
quadrupole pattern of the form uq = VqT ■ f, where Vg is a 
traceless, symmetric 3x3 tensor. They have attributed this 
"VELMOD quadrupole" residuals to uncertainties in modeling 
the IRAS density field beyond 3 000 km s" 1 . In particular they 
found that the most important sources of uncertainties were the 
smoothing procedure (based on a Wiener Filtering technique) 
and the shot noise errors, with the errors deriving from having 
ignored mass inhomogeneities beyond 12,000 kms" 1 playing 
a minor role. Using PSCz instead of 1.2-Jy catalog reduces the 
shot noise errors in the density field beyond 3000 kms" 1 and 
thus should improve the agreement with the Mark III velocities. 
In other words one would expect the quadrupole of the resid- 
ual velocity field Mq 6 - uf 2 , where one only takes in account 
the mass beyond 3000 km s" 1 , to be similar to the "VELMOD 
quadrupole" of WSDK. 

In Fig. (15) we compare the VELMOD quadrupole at s = 
2000 km s" 1 (rescaled from j3 = 0.492 to (3 = 0.6, upper panel) to 
the Unf.-uf 2 quadrupole residual (lower panel). The VELMOD 



-*0.6 "1.2 1 

quadrupole reaches its maximum value at (f « 165°, b s 



55°) 



and in the opposite direction of the sky. The Mq 6 - u\ 2 resid- 
ual quadrupole exhibits a rather different pattern, reaching its 
maximum amplitude at (I w 156°, b k, -74°). The value of the 
quadrupole components of 1.2-Jy and PSCz as well as those of 
the VELMOD and UQ b -uf 2 quadrupole residuals computed at 
a distance of 2000 kms" 1 are shown in Tab. (4). This com- 
parison is at best qualitative since (f) VELMOD quadrupole is 
computed in real space while both IRAS quadrupoles refer to 
redshift space insted, and (if) the smoothing scheme used in 
WDKS is different from with ours. 

7. CONCLUSIONS. 

We have compared the model velocity and density fields of 
the PSCz survey to those derived from the sparser 1 .2-Jy sam- 



ple in redshift space. Our model predictions are based on 
the Nusser & Davis (1994) spherical harmonics decomposition 
technique and assumes Zeld'ovich approximation and linear bi- 
asing. 

The IRAS overdensity fields are reconstructed with a smooth- 
ing proportional to the 1.2-Jy inter-galaxy separation and ex- 
hibit the same general features corresponding to the known cos- 
mic structures in the local universe. The model velocity field, 
computed in the LG frame, has a dipole-like appearance, with 
galaxies moving away from us toward Perseus-Pisces region 
and approaching the LG from the opposite direction. 

The analysis of the velocity multipoles has revealed a signif- 
icant discrepancy between the PSCz and 1 .2-Jy velocity mono- 
poles in the range 1 800 < s < 4000km s" 1 in which PSCz radial 
velocities are systematically larger than 1.2-Jy ones, that cannot 
be explained by uncertainties in the IRAS selection functions 
or incompleteness of the PSCz catalog at faint fluxes. Also, 
the mismatch cannot be ascribed to errors derived from shot- 
noise and treatment of unobserved areas since both have been 
accounted for with the help of realistic mock IRAS catalogs and 
the extensive application of bootstrap resampling techniques to 
both IRAS datasets. 

Both IRAS dipole and quadrupole moments are in good 
agreement within the sampled volume apart from a small mis- 
match in the u\o, M2-2 and M22 quadrupole components beyond 
s = 7 000 kms" 1 . The agreement between the dipole moments 
and their consistency with the velocity dipole computed by BD- 
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FIG. 15. — VELMOD (upper panel) and »o 6 _ »f2 ower panel) quadrupole 
residuals at s = 2 000 km s . Stars denote outflowing from us (LG frame); cir- 
cles denote infall to us. The size of the symbol is proportional to the amplitude 
of the flow (notice the different velocity scales in the two panels). The max- 
imum amplitude of the "VELMOD quadrupole" (top panel) is 181 kms -1 , 
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SLS from the ORS catalog suggests that deviations from ho- 
mogeneous linear biasing prescriptions have to be small and 
cannot be advocated to explain the aforementioned mismatch 
between the IRAS monopoles. 

Another remarkable consequence of the agreement among 
the various dipoles is that the 1 .2-iy-Marklll residuals cannot 
be alleviated by using the PSCz catalog instead of the 1 .2-Jy or 
the ORS ones. In fact, this suggests that the origin of the 1.2- 
iy-Mark III dipole resides in the Mark III calibration procedure; 
a suggestion corroborated by the fact that the dipole mismatch 
disappears when adopting the alternative VELMOD procedure 
to calibrate -Marklll velocities (WDSK). The VELMOD anal- 
ysis requires a sizable external quadrupole component within 
3000 kms" 1 . This quadrupole accounts for uncertainties er- 
rors in the 1 .2-Jy model density field beyond that distance. We 
found that using the PSCz sample provides part of the required 
external tidal field. 

As a final remark, we would like to stress that the present 
analysis illustrates once more the importance of spherical har- 
monics decomposition techniques to compare observed and 
model density and velocity fields obtained from all-sky sam- 
ples, to reveal possible data-model or model-model inconsisten- 
cies and to assess model adequacy. These techniques will prove 
very useful when the next generation of all-sky redshift surveys 
and peculiar velocity catalogs [e.g. the 6dF (http:www.mso.anu. 
edu.au/6dFGS) and 2MRS (http:cfa-www.harvard.edu/huchra/ 
2mass) datasets] will become available. 
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Table 4 

Quadrupole Components at 2 000 km s~' 



VELMOD PSCz 1.2-Jy PSCz - 1.2-Jy 
[kms -1 ] [kms -1 ] [kms -1 ] [kms -1 ] 



Vq (x,x)... 


45 


57 


69 


-13 


Vq (y,y)... 


44 


-22 


-15 


-7 


Vq (z,z)... 


-89 


-35 


-54 


20 


Vq (x,y)... 


18 


-150 


-156 


6 


Vq (x,z)... 


138 


57 


48 


10 


Vq (y,z)... 


-29 


-106 


-102 


-5 
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NOTE.-(l) We have rescaled the "VELMOD quadrupole" compo- 
nents presented in WSDK's Tab. (2) to the numerical value of (3 = 0.6. 



